Why some Iron-based superconductors are nodal while others are nodeless 
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The symmetry of the order parameter in iron-based superconductors, especially the presence or 
absence of nodes, is still a question of debate. While contradictory experiments can be explained 
by appropriately tuned theories of nodeless superconductivity in the iron-arsenide compounds, for 
LaOFeP all experiments clearly point to a nodal order parameter. We put forward a scenario 
that naturally explains the difference between the order-parameter character in these two sets of 
compounds, and use functional renormalization group (fRG) techniques to analyze it in detail. Our 
results show that, due to the orbital content of the electron and hole bands, nodal superconductivity 
on the electron pockets (hole pocket gaps are always nodeless) can naturally appear when the third 
hole pocket which lies at wavevector (n, n) in the unfolded Brillouin zone is absent, as is the case 
in LaOFeP. When present, the third hole pocket has overwhelming d xy orbital character, and the 
intra-orbital interaction with the d xy dominated part of the electron Fermi surface is enough to drive 
the superconductivity nodeless (of s form). However, in its absence, pair hopping, inter-orbital, 
and electron-electron intra-orbital interactions render the gap on the electron pockets softly nodal. 

PACS numbers: 74.20.Mn, 74.20.Rp, 74.25.Jb, 74.72. Jb 



After two years of intense research in the physics of 
the new iron- based superconductors [l[ , the symmetry of 
the order parameter is still far from settled. Theoret- 
ically, the current opinion converged on a s± nodeless 
order parameter that changes sign between the electron 
(e) and hole (h) pockets, which comes out of both the 
strong and the weak-coupling pictures of the iron-based 
superconductors However ^most of these theories 

are either phenomenological 0, Hf , or use models that 
take into account only a part of the Fe orbitals present 
in the material. The most quantitative approaches give 
anisotropic gaps around the electron Fermi surface which, 
at their smallest value, are close but do not cross zero Q . 



The experimental situation is more controversial. In 
the 122 compounds, ARPES points to the existence of 
nodeless isotropic gaps on the hole Fermi sur faces [Hill; 
on the electron Fermi surface, on which ARPES also 
shows an isotropic gap, the data is less trustful. NMR 



spin relaxation time shows no coherence peak [131 ] and 
exhibits a T 3 power-law right below the transition tem- 
perature up to 0.1 — 0.2T C , reminiscent of the cupratcs 
nodal behavior. This can be explained by fine-tuning 
an s± order parameter [lij ]. Penetration depth data on 
the Fe-As compounds display power-law |l5| and expo- 
nential [l6j behavior. Thermal conductivity shows no 
residual k/T intercept, although its in- field dependence 
points to a large gap anisotropy [13 EH- By contrast, for 
LaOFeP, the penetration depth is extremely linear with 
temperature, and the thermal conductivity shows a large 
residual intercept, both strong characteristics of nodal su- 
perconductivity 19|, 2(J . This is even more puzzling since 
both materials have similar Fermi surface topology plj ]. 
In this article, we offer an explanation for the puzzling 




FIG. 1: 2D Band structure for LaOFeAs (a) and LaOFeP 
(b). (Inset: Brillouin Zone). The tight-binding model for 
LaOFeAs is contains in-plain hoppings up to fifth nearest- 
neighbors [jj; the parameters are varied for LaOFeP accord- 
ing to the different pnictogen height parameters given in [IoT | . 
The dashed horizontal lines denote the Fermi Level for the un- 
doped compound. The electronic structure looks very similar 
in both systems. The major difference is the d X 2-Y 2 dom- 
inated band crossing the Fermi level in (a), but not in (b). 



Being still away from the Fermi level, the 3z 
band is shifted up in (b) compared to (a). 



dominated 



difference between the order parameter character in the 
As and P-based compounds. Using functional renormal- 
ization group (fRG) on a 5-band orbital model of the 
iron-based superconductors with orbital interactions, we 
find that the gap on the e-pockets can undergo a nodal 
transition if the third hole pocket at (jr, 7r) is absent. 
Both ARPES and numerical simulation data show this 
to be true in the case of LaOFeP [22], [23| . By using an 
extended orbital model due to Kuroki et al. [H], [Toj , we 
find that the third h-pocket is overwhelmingly composed 
of d xy orbital (or dx 2 -Y 2 m a 45° rotated basis along 
the Fe-As bonds). When present, the third h-pocket at 
the (tt, 7r) point in the unfolded Brillouin zone (Fig. 
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FIG. 2: 5 pocket scenario for As-based compounds, a) Plot of the superconducting form factor gap versus the patching indices 
(momenta) shown in b). The gap on the outer h-pocket at F is smaller than of the inner h-pocket and of the same order as 
the M pocket gap. The gap on the e-pockets is very anisotropic but nodeless and of opposite sign from h-pocket gap. cl)-c3) 
Orbital weight distribution on the different pockets (not shown is the outer h-pocket at V, which is similar to c2) shifted by 
90 degrees assuring orthogonality of the band vectors). Dashed lines indicate most relevant scattering contributions for the 
dominating Ui intra-orbital interaction, d) Leading pairing instability eigenvalues cf^ a (,)(A C ) from the orbital decomposition 
matrix (clca) (cbCb) in the Cooper channel. Most relevant weights are on the diagonal and off-diagonal contributions from 
the dxz,YZ orbitals as well as the diagonal contribution of d X 2_ Y 2 - e ) Flow of leading instability eigenvalues (charge density 
wave (CDW), Pomeranchuk instability (PI), SDW, and SC). SDW fluctuations are highly relevant, leading (s ± ); sub-leading 
(d-wave) SC instabilities as well as SDW diverge in close proximity to each other, where the divergence scale is of A c w 0.08eV . 



scatters through intra-orbital interactions strongly with 
the d xy part of the e-pockets around the X point to form 
an phase. This enhances the already present su- 
perconductivity between the T point h-pockets and the 
e-pockets. The gap on the e-pockets is found to be very 
anisotropic. Upon changing to the P-based structure, 
the absence of the third h-pocket allows for pair hop- 
ping, inter-orbital interactions, and electron-electron (e- 
c) scattering to drive the already anisotropic gap on the 
c-pockcts nodal. We use a two-dimensional tight-binding 
model developed by Kuroki et al. [B| to describe the band 
structure of the 1111-typc iron-based superconductors: 



H 



EE 

k.s a, 6—1 



,(k)c k 



(1) 



Here c's denote electron annihilation operators, and a, b 
the 5 d Fe orbitals and s the spin indices. While the main 
electronic structure of P-based and As-based compounds 
is very similar, there are certain important differences. 
Fig.Q]shows the band structure of LaOFeAs and LaOFeP, 
where the latter is obtained by adjusting the parameters 
in Q according to the changed pnictogen height from As 
to P [13]. In the vicinity of the Fermi surface, the most 
notable difference is the presence or absence of a broad 
dx^-Y 2 (dzj^-dominated band at M = (tt,tt), in agree- 



ment with ARPES data. To account for this difference, 
we use a 5 pocket scenario for the As-based and a 4 pocket 
scenario for the P-based compounds. We choose to com- 
pare and analyze two generic cases, with or without the 
h-pocket at (7r,7r) as corresponding to the As-based 1111 
(122) and the P-based compounds, respectively. 
The interactions in the orbital model are given by: 



H in t = 



^E 



a<b s,s' 



'h aii' c Lt C iai C i b± C ib J [ ) 



where ni^ as denote density operators at site i of spin 
s in orbital a. We consider intra- and inter-orbital 
interactions Ui and U2 as well as Hund's coupling Jh 
and pair hopping J pa ir- In what follows, we choose 
a physical interaction setting dominated by intra- 
orbital coupling, U\ > U-2 > Jh ~ Jpair, and choose 
U x = 3.5eV,[/ 2 = 2.0eV,Jn = J pair = 0.7eV @. (Even 
though the interaction scales are chosen relatively high, 
the bare effective interaction scale, taking into account 
different orbital weights, does not exceed ~ 2eV, versus a 
kinetic bandwidth of ~ 5eV.) This is a rather simplified 
picture as these scales would in general also depend on 
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FIG. 3: 4 pocket scenario for LaOFeP. a) Leading s superconducting form factor gap, with Fermi surface patches given in b); 
the d X 2-Y 2 -dominated h-pocket at M is absent. The h-pockets at V are gapped and rather isotropic, with a smaller gap on the 
outer h-pocket. The e-pockets show a huge anisotropy, being nodal on the inner pockets tips indicated by dashed arrows, c) 
Important processes leading to nodal channels. Being subleading in the presence of the M hole pocket, J pa ir scattering processes 
from the dxz,YZ dominated V pockets to the e-pockets and electron-electron intra-orbital processes between e-pockets now 
contribute to nodal propensity, with the orbital weight distributions shown in fl)-f2). d) The orbital decomposition matrix of 
the SC instability show increasing weight on the dxz,YZ orbitals with a less relevant diagonal d X 2_ Y 2 contribution, e) Flow of 
leading instability eigenvalues (notation as in Fig. [2^); nodal a is favored; the SDW fluctuations are comparably weak. The 
divergence scale A c w 0.0065el/ is much less than for the As-based materials. 



the different orbitals. However, for the two scenarios rep- 
resentative for As- and P-based compounds, we checked 
that the main features are stable under variation of these 
parameters. From the band structure point of view, it 
should also be noted that the d^ z 2_ r i -dominated band 
moves towards the Fermi level for the P-based compound 
(Fig. [1]). However, this band only plays a marginal role 
since no other relevant band has this orbital content 
and hence any scattering to other bands would be 
governed by sub-leading inter-orbital interactions. Wc 
defer a refined discussion to a later stage 01. Using 
functional renormalization group (fRG) fed - Eil ] , we 
study how the renormalized interaction described by 
the 4-point function (4PF) evolves under integrating 
high energy fermionic modes; the flow parameter is 
the IR cutoff A approaching the Fermi surface, i.e. 
y A (fc 1; n l5 fc 2 , n 2 ; fc 3 , n 3 ; fc 4 , n 4 )4 4n4S 4 3 „ 3? c fc2ii2S c fcitii -, 
with ki to &4 the incoming and outgoing momenta. 
Due to the spin rotational invariancc of interactions, we 
constrain ourselves to the S z = subspacc of incoming 
momenta k\,k2 (and outgoing ^3,^4) and generate 
the singlet and triplet channel by symmetrization and 
antisymmetrization of the 4PF Va [27| . The starting 
conditions are given by the kinetic bandwidth serving as 
an UV cutoff, with the bare initial interactions for the 
4PF. The diverging channels of the 4PF under the flow 
to the Fermi surface signal the nature of the instability, 
and A c serves as an upper bound for the transition 



temperature T c . For a given instability characterized by 
some order parameter Ok (the most important example 
of which is the SC instability (5| c = in our 

case), the 4PF in the particular ordering channel can 
be written as Z)fe, p VA.(fe,p)[6^6 p ] (29|. Accordingly, 
the 4PF in the Cooper channel can be decomposed into 
different eigenmode contributions 

F A sc (fc,-fc,p) - £ Cl sc (A)/ sc W/ SC '>), (3) 

i 

where i is a symmetry decomposition index, and the lead- 
ing instability of that channel corresponds to an eigen- 
value cf C (A) first diverging under the flow of A. / SC ' z (fc) 
is the SC form factor of pairing mode i which tells us 
about the SC pairing symmetry and hence gap structure 
associated with it. In fRG, from the final Cooper chan- 
nel 4PFs, this quantity is computed along the discretized 
Fermi surfaces, and the leading instability form factors 
are plotted in Fig. [^i and [3^. 

As-based compounds: For the As-based setting, we find 
that the instability is the leading instability of the 
model at moderate doping. The setup resembles the sit- 
uation studied in Q, which, as an additional check, we 
also studied with a more detailed tight-binding structure 
beyond 5th next-nearest neighbors. We likewise find a 
nodclcss s^ SC leading instability. However, we can iden- 
tify the M h-pocket to play a major role in contributing 
to the SDW fluctuations and to support the full gapping 



4 



of the s-wave as well as sign-change from hole to elec- 
tron pockets (Fig. [2]). In particular, we study the orbital 
content in detail and analyze how the pairing instability 
distributes over the different orbitals (Fig.[2Ji). For this, 
we consider the 4PF in orbital space, 



v: 



orb 

c,eZ— >a,b 



^ | VA.(feij n\\ k 2 ,n 2 ; k 3l n 3 ; k 4 ,n 4 ) 



n± , . . . ,714 — 1 



X < ni ( k l)u*bn 2 {k 2 )u cn;t (fc 3 ) (**)}, (4) 

where the it's denote the different orbital components of 
the band vectors. The matrix shown in Fig. [3Ji gives 
the leading eigenvalue contributions of V^ Q b ^(k,p) = 

( C I a c -k a)( c p,b c -p.b)i m the Cooper channel of (jl]). As 
in we decompose it into different form factor contri- 
butions Ei^ a , 6) (A)/g(fe)*/gS(p), where the lead- 
ing eigenvalues at A c for different (a, b) are plotted in 
Fig. [2ji and[3ji. We observe a leading contribution in the 
diagonal part of the dx 2 -Y 2 orbital, which is strongly 
linked to the scattering contributions from the M h- 
pockct to thc dx 2 -Y 2- dominatcd parts of thc e-pockcts 
(Fig. (2). The main scattering processes are intra-orbital 
scattering between the d xz (or d yz ) orbitals-dominatcd 
parts of the e- and T h-pockets. These processes favor 
an SC instability, as was found in Q. However, there 
is another process in the case of As-based superconduc- 
tors: Along the r «-> X path, the e-pocket has a high 
concentration of the d X 2 -Y 2 orbital. This part of the 
c-pocket then scatters strongly with the 3rd h-pocket at 
the M-point, which is entirely made of the dx 2 -Y 2 or- 
bital band. The intra-orbital repulsion on them prefers 
an s ± -type pairing between the M h-pocket and the e- 
pockct, which reinforces the already present between 
the r h-pockets and the X (X 1 ) c-pockets. With the 
commonly-used assumption that t/i n t ra is the dominant 
interaction, it then seems likely that the 3 h-pockets dis- 
play a gap of identical sign: the two T-pockets, which are 
not nested with each other and hence can have the same 
gap sign and are of different orbital content than the M 
h-pockct. This then allows all these gaps to have the 
same sign. However, the e-pockets contain contributions 
from all 3 mainly relevant d orbitals, and therefore scat- 
ter strongly through Z7i n t ra with all 3 h-pockets, which 
enhances the character of the gap. The presence of 
the 3rd hole pocket is also responsible for the strong SDW 
signal, being due to the fact that the nesting wave- vector 
M o X is the same asTnl. Below we will indeed see 
below that the absence of this pocket weakens the SDW. 

P-based compounds: Moving to the P-based compound 
case changes the physical picture even qualitatively. As 
shown in Fig. [3j we find a nodal scenario for the 
P-based compounds, with lower critical divergence scale 
A c ~ T c and less SDW-type fluctuations. The absence of 
SD W order in the P-based compounds is a well-known ex- 
perimental fact. Here, the absence of the M hole pocket 



removes the intra-orbital scattering contribution to the 
electron pockets. This gives way to previously sublcad- 
ing scattering channels like e-e scattering between the 
dx 2 -y 2_ dominatcd parts of the e-pockets and pair hop- 
ping from the h-pockets at T to the e-pockets. They tend 
to drive this part of the e-pockets nodal (Fig. [5]:). As an- 
other check, we applied constant band interactions to the 
As-based and P-based setups to analyze the pure Fermi 
surface topology unaffected by orbital-specific interaction 
effects. This means that we do not take into accout the 
orbital content at a Fermi surface point, but choose a 
general set of couplings g that is only specified by the 
type of intra- or inter-pocket scattering 25, 26|. We find 
that the intensity of (0, 7r), (tt, 0) SDW fluctuations is in- 
deed strongly influenced by the scattering contributions 
M -B- X. Constant band interactions show very isotropic 
gaps on the e-pockets [H, [26J. As a consequence, the 
electron gap anisotropy, rendering the electron gap node- 
less or nodal, is a phenomenon of orbital interactions, in 
accordance with above elaborations. Furthermore, we 
varied the intra-orbital interaction scale on the dx 2 -Y 2 
orbital and the J pa ir scale in the P-based scenario. In- 
creasing these scales, it triggers the propensity to form 
nodes on the electron pockets. Furthermore, under this 
variation, the electron pocket part of the form factor can 
be considerably shifted to being more or less nodal. 

These tendencies all appear to be consistent with ex- 
periment: we find (i) a lower divergence scale and hence 
lower critical temperature, (ii) significantly enhanced low 
energy density of states in the superconducting phase, 
and (iii) reduced SDW type fluctuations, which, even at 
pronounced nesting, are insufficient to drive the system 
to a leading magnetic SDW instability [l^-tlH. The ab- 
sence of the hole pocket at M also manifests itself in the 
orbital decomposition of thc pairing instability Fig. [3J1: 
The diagonal contribution of d X 2 ~Y 2 -i m comparison to 
the dxz,YZ, is significantly reduced. 

We find that the relevance of the broad band at the 
unfolded M point plays the major role to explain the 
drastic change of properties from the As-based to the 
P-based 1111 compounds, rendering thc former nodclcss 
and the latter nodal. The nodes that appear in the P- 
based compounds are driven by anisotropy and are not 
given by any new SC pairing symmetry, which remains 
of s^ type. A more detailed work on the technical back- 
ground and broader scope of this work is in progress . 
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